Screened electrostatic interactions between clay platelets. 
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An effective pair potential for systems of uniformly charged lamellar col- 
loids in the presence of an electrolytic solution of microscopic co- and counteri- 
ons is derived. The charge distribution on the discs is expressed as a collection 
of multipole moments, and the tensors which determine the interactions be- 
tween these multipoles are derived from a screened Coulomb potential. Unlike 
previous theoretical studies of such systems, the interaction energy may now 
be expressed for discs at arbitrary mutual orientation. The potential is shown 
to be exactly equivalent to the use of linearized Poisson-Boltzmann theory. 

I. INTRODUCTION 

While the mesostructure, stability and phase behaviour of charge-stabilized dispersions 
of spherical colloidal particles are by now reasonably well understood, both experimentally 
and theoretically |]-|| , the picture is much less clear in the case of lamellar colloids, of which 
clay dispersions are a pre-eminent example Q . This is partly due to the high degree of poly- 
dispersity, the irregular shapes, and the extreme anisotropy of the thin lamellar particles of 
naturally occurring clay suspensions. Such complications render an unambiguous interpreta- 
tion of experimental data, eg. from small angle X-ray or neutron diffraction measurements, 
very difficult, while posing a practically insurmountable challenge to the theoretician at- 
tempting a Statistical Mechanics description. Even for the widely studied synthetic model 
system of Laponite, made up of nearly monodisperse, disc-shaped platelets, there is no 
consensus among experimentalists as to the structure, gelling behaviour and rheology of 
semi-dilute suspensions ||-||, while attempts at a theoretical description, or simulations of 
this model system are in their infancy [ n^fnj . The main reason for the latter state of affairs 



is that a realistic model for the effective interaction between a pair of arbitrarily oriented 
charged circular platelets, generalizing the isotropic DLVO potential |IJ between spherical 
colloids, is not available. Only in the simplest case of two coaxial, uniformly charged discs, 



has a screened Coulomb interaction been worked out within linear |E| and non-linear fl3 
Poisson-Boltzmann (PB) theory. 

The Molecular Dynamics (MD) simulations of ref. [11] were based on a site-site inter- 
action model, generalizing the "Yukawa segment" representation used earlier to simulate 



suspensions of charged rods [14]. Within such a representation the charge distribution on 
a rod or a platelet is discretized into v interaction sites, each carrying a fraction 1/v of 
the total charge; sites on different particles interact via a screened Coulomb potential. This 
Yukawa segment model is very computationally intensive, since the total interaction between 
two particles involves v 2 contributions. It also carries a degree of arbitrariness in the choice 
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of the number of sites v, which for practical reasons must be taken to be much smaller than 
the number of elementary surface charges (typically 10 3 for Laponite) carried by individual 
platelets. 

This paper examines the continuous version of the Yukawa segment model, and derives a 
multipolar expansion of the effective, screened Coulomb interaction between two uniformly 
charged platelets of arbitrary relative orientations. The resulting anisotropic effective pair 
potential is shown to be accurate for centre-to-centre distances larger than the radius of 
the platelets, and should hence provide a useful tool for theoretical investigations of the 
structure and sol-gel transition in semi-dilute clay dispersions. 



II. THE YUKAWA SEGMENT MODEL 

Consider a suspension of N p infinitely thin circular platelets per unit volume, of radius 
a, and carrying a uniform surface charge density a = Ze/na 2 , where Ze (< 0) is the total 
charge on a platelet. The platelets, together with microscopic co- and counterions, are 
suspended in water. Since the present study focuses on mesoscopic lengthscales, of the 
order of a ( typically a ~ 15nm for Laponite), one may neglect the molecular nature of 
water which will be regarded as a continuum of dielectric constant e. Within linearized PB 
theory, the effective interactions between platelets are always pairwise additive |1| , and the 



screening of electrostatic interactions by the co- and counterions is uniquely characterized 
by the Debye screening length: 

where the sum is over all microion species, n a and z a are the concentration (number density) 
and valence of ions of species a, and €q is the permittivity of free space. Building on the 
linearity of the theory, the Yukawa segment model assumes that each infinitesimal area ds 
(or "segment" ) on a uniformly charged disc generates a screened Coulomb potential: 

-e Kr . (2) 
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The corresponding pair potential between two infinitesimal areas ds (around r) and ds 
(around r ) on two discs is then 

„(| r _ r '|) = a2dsds ' - e-^'l, (3) 
Ane e\r - r | 

and the total pair interaction between the discs is obtained by integrating (|3|) over the 
surfaces of the two discs. However, for arbitrary orientations of the discs, this leads to 
intractable expressions involving multiple integrals. 

Instead, by analogy with electrostatic interactions between extended charge distributions, 
a systematic multipolar expansion of the screened Coulombic interaction will be sought. The 
derivation of such an expansion requires some care, because the basic screened Coulomb 
potential (|2|) does not satisfy Poisson's equation, except in the bare Coulomb limit, where 

K — ► 0. 
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III. THE POTENTIAL AROUND A SINGLE PLATE 



To get a feeling for a multipolar expansion involving screened, rather than bare, Coulomb 
interactions, consider first the potential due to a uniformly charged disc, along the axis of 
the disc. Using cylindrical coordinates, with the z— coordinate along the axis of the disc (cf. 
Fig. [I]), the screened potential along that axis (ie. for a radial coordinate p = 0) is simply 

V(p = 0,z) = / RdR 

KH ' 47re eio ^z 2 + R 2 



= ^[[1_ e-^V^tT- 1 ]] (4) 

which is easily expanded in powers of (a / z) according to 

7e 00 / 1 \ n+1 

9(p = 0,z) = j e~^\ Y,[An + K n }[-\ , (5) 

4vre e n=0 \\z\J 

where the coefficients A n and K n are listed in Table [l|. 

Several points are to be noted about this expansion. First, the corresponding expansion 
for the bare Coulomb potential is correctly retrieved by taking the limit k — > 0; this amounts 
to setting all K n = in eq.(|5|), leaving only odd powers of 1/z in the expansion, since all odd 
coefficients A n vanish. This is an obvious consequence of the space reflection symmetry of 
the uniform charge distribution on a disc, which implies that only even multipole moments 
exist. However, for the screened Coulomb potential, terms with even powers of 1/z appear in 
the expansion, which would correspond to odd multipoles (dipole etc.) in the bare Coulomb 
case. 

The second remark is that the expansion @ also follows from the exact potential due to 
a uniformly charged disc, and its associated electric double-layer of co- and counter ions, as 
calculated within linearized PB theory ||16|| , namely: 

2Ze r°° „ s „ p-W^ 3 - 



*(p,z) = / J 1 (ka)J (kp)-===-dk. (6) 

ae e Jo \Jk 2 + k 2 

This agreement between the expansion in eq.(||D and the expansion of eq.(g) for p = is an 
illustration of the exact equivalence between the Yukawa segment model, and a full linearized 
PB calculation of the effective potential generated by a charged particle of any shape and 
its associated electric double layer. 

The final point concerns the generalization of the expansion (|J) away from the z— axis, 
ie for p 7^ 0. In the bare Coulomb case (k = 0), the coefficients A n in table 1 may be 
immediately carried over to spherical polar coordinates (r, 9, <p) to write down an expansion 
of the potential due to a uniformly charged disc in even Legendre polynomials (the potential 
is independent of the azimuthal angle <p): 



Ze 00 1 
*M) = 4^ £ A^Pn(cos0). (7) 

t''tO c n =0; even ' 



However, the presence of the exponential screening factor in the expansion (|5|) prevents a 
similar straightforward generalization to off-axis conditions in the screened Coulomb case 
(k ^ 0). For this reason the multipolar expansion must be re-examined more carefully for 
the Yukawa segment model. 
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IV. SCREENED MULTIPOLAR EXPANSION 



The multipolar expansion of the total potential $f(r) due to a uniformly charged disc, 
with each infinitesimal surface element generating the screened potential (2) (Yukawa seg- 
ment model), may be derived along lines nearly identical to the classic calculation for un- 
screened charge distributions |T7|]. Clearly 

*(r) = f <j)(\r-s\)ds (8) 
J s 

where the integral is over the surface S of the disc. The potential </> may now be expanded 
as a Taylor series about the centre of the disc (s = 0), ie.: 

*(r) =<r f ds\ <P{r) - J2 s a V a <P{r) + ^ E s a s p V a V p<j>{r) + ■■•!, (9) 

JS { l - a/3 J 

where the sums are over all three Cartesian coordinates of the vector s G S. All odd terms 
(eg. the dipolar term) vanish by symmetry. This leaves only the even terms: 

*(r)=^(r)+* (r)+**(r) + ..., (10) 

involving the total charge Ze of the disc, its quadrupole tensor Q, its hexadecapole tensor 
J>, etc.: 

*z(r) = ZeT K (r) 

*Q(r) = ^QapT^ir) (11) 

where the Einstein convention of summation over repeated indices has been adopted. The 
tensors T K are: 

7V. = W7 ...(— — ) (12) 

while the Q a /3, QapyS are the Cartesian components of the 2 nd rank quadrupolar and 4 th rank 
hexadecapolar tensors. For a uniformly charged disc the quadrupolar tensor is given (in a 
frame where the z-coordinate is along the axis of the disc) by 

Q = — ssds 
= e Js 

( -Q \ 

= -Q (13) 

V 0/ 

where Q = —Za 2 /4. Note that contrary to the bare Coulomb case, Q cannot be chosen 
traceless in the screened case, because the tensors T£a are themselves not traceless. The 
Cartesian components of the 4 th rank hexadecapole moment are defined by: 
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$a/3 7 <5 = - / SaSpS-fSsds. 

e Js 



(14) 



For a disc, choosing the z— coordinate along its axis, the only non-zero components aj3 r )8 
are the two diagonal components $ xxxx and and those in which x and y both appear 

twice. Explicitly 



^xxxx ^yyyy 



Za 4 



^xxyy ^xyxy 



Za 4 _ $ 
~24~ = 3' 



(15) 



The calculation of the tensors is considerably lengthier for the screened than for the 

bare Coulomb interaction. Some details are given in Appendix A. 

In spherical coordinates, the total potential due to the uniformly charged disc, up to 
hexadecapolar order, may finally be written as: 



47reoe 



Z + Q K 2 (cos 2 - 1) + — k 4 sin 4 



24' 



+- 



1 

1 

+ ^r 



Qk(3 cos 2 9 - 1) - -k 3 ( sin 2 - — sin 4 
2 v ; 3 V 32 

5(3 cos 2 - 1) + -k 2 f 1 - 6 sin 2 9 + — sin 4 9 
2 3 V 8 

-k(35 cos 4 0-30 cos 2 + 3) 



-(35 cos 4 0-30 cos 2 + 3) 



(16) 



Returning to cylindrical coordinates, this expression reduces on the z— axis to: 



*(p=0,z) 



< 

47re e 




(17) 



which coincides with the expansion ([5]), with coefficients given in Table 1, up to order a 4 (the 
higher powers of a in the coefficients of 1/z 4 and 1/z 5 corresponding to higher order multipole 
moments). To illustrate the convergence of the multipolar expansion of the potential to 
hexadecapolar order, in figs. ^| and |3| the potential given by flTE| ) is compared with the 
linearized PB potential of eq.(^) and an explicit numerical integration over a discretized 
charge distribution (with v = 7841 sites), both along the z— axis and in the xy— plane. 
The numerical integration was carried out to check its accuracy against the exact potential 
(§), since it will be the only available test of the multipolar expansion of the screened pair 
interaction between two platelets at arbitrary orientations, to be discussed in section VI. 
Agreement is seen to be excellent, down to around z(p) ~ a, where the multipolar expansion 
diverges dramatically, as higher order terms in (l/r) n start to dominate. 
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V. EFFECTIVE INTERACTIONS BETWEEN TWO PLATES 



The results of the previous section for the effective potential around a single platelet may 
now be used to determine the potential energy of interaction (or effective pair potential) of 
a pair of discs, at arbitrary relative orientations, as shown in Fig.|]. The interaction energy 
Vab is formally expressed by integrating the screened electrostatic potential arising from 
one disc (A) defined by eqs.fllQf) and (0), over the surface charge distribution of the second 
disc (£): 



Vab^Oa^bAaAb) = I crds^ A (r + s) 

J Sb 



(18) 



The potential ty A (r + s) may be expanded in a Taylor series in powers of s, about the centre 
of disc B (s=0), along the lines of eq.(pf), with the result: 



V AB = aj s ds !^ A (r) + sf V„* A (r) + ±tg s B p V a V^ A (r) + ■ 



(19) 



Now inserting the expansions of the potential \& (r) (10), and its derivatives, into (19) the 
electrostatic pair potential is given by 



Vab = ea ds 1 
Js B 



rprz ry A i ^ rpn f^iA i ^~ rpn (f\A i 
rpK yA I rpK f^A _j_ ^~ rpK f^A i 



t~T~tK fyA I Tif; _L 
J «/3 Z; i* 27 J a/375 V 7 <5 t " 



+ ■ 



(20) 



Integrating over the surface, Sb, of disc B naturally introduces the multipole moments of the 
surface charge distribution of that disc into the expression ( [2~0"|) for Vab, which may hence 
be cast in the form 



v AB (r, & Ai e B , <p A , <Pb) = vjg + tyfg + v%) + vfif + (vjg + v%) + 



rQZ^ 



(21) 



where 



Vf£ = e 2 Z A T K Z A 



V 



ZQ 



z a t: q 



AB 



V AB 



2! 



Ma/3 



2!2! 



Qa(3^af3ysQj8 



V 



AB 



yArriK fcB 



(22) 
(23) 

(24) 

(25) 



This expansion is consistent up to order 1/r 5 with the corresponding series for pure Coulom- 
bic interactions (k, = 0). The rather cumbersome expressions for the screened interaction 
tensors are given in Appendix A, while details of the summations over Cartesian coordinates, 
implicit in eqs.(|23]-E5l), are described in appendix B. 
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VI. RESULTS FOR SPECIFIC GEOMETRIES 



The detailed behaviour of the effective pair potential Vab defined by eqs. , will 



now be examined, as a function of the centre-to-centre distance between, and the mutual 
orientations of, the two discs. The relevant variables are 

• r, the separation of the centres of masses of the two discs 

• {@a, @b} the polar angles of the two discs 

• A0, the difference in the azimuthal angles of the 2 discs 

• k, the inverse screening length determined by the co- and counterions, 

Rather than specialize to the physical parameters particular to Laponite (Z ~ —1000, a ~ 15 
nm), in the figures which follow all distances shall be scaled by the disc radius a, which 
provides a convenient lengthscale, and all energies by the bare Coulombic energy of two 
charges Ze a distance a apart, ie. by Z 2 e 2 /47ieoea. Thus these results apply to the most 
general uniformly charged disc. 

The case of two coaxial plates is illustrated in Fig. |5|, where the various contributions 
to the energy Vab are plotted as a function of the distance between the two plates. The 
exact result in this simple geometry is known within LPB theory [I2|], and used to test 



the convergence of the multipolar series. As expected, the charge-charge contribution ( |22| ) 
dominates for large spacings, but the contributions of the higher order terms become very 
significant for spacings less than the disc diameter, 2a. 

In Fig. P the energy is plotted as a function of the separation between two discs lying in 
the same plane, like two coins on a table. In this geometry all multipolar contributions are 
positive (repulsive), and the enhancement of the total energy over the charge-charge term 
is more pronounced. Included in this figure is the potential energy obtained by integrating 
over the Yukawa potential over two discretized site-charge distributions, where each disc has 
v = 7841 point charges, each carrying charge Zejv. The agreement of this calculation with 
the multipolar expansion is excellent, down to the distance of closest approach r = 2a below 
which point the discs intersect. 

In Fig. [7| the energy is plotted as a function of the separation of two discs in a T-shaped 
configuration (9a = 0, 9b = tt/2), a geometry favoured by the quadrupoles, as also observed 
in ref. ]I(J where a purely quadrupolar disc model was studied. The agreement with the 



discretized Yukawa segment calculation is again excellent, except at very close centre-to- 
centre separations, where the multipolar expansion is expected to collapse. 

In the next plot, Fig|| the behaviour of the pair potential is examined as one disc 
slides over the other at constant altitude (z = 1.5a in this plot), with both discs parallel 
to each other. As expected the energy goes through a maximum at the distance of closest 
approach, but detailed structure is seen in the total energy, as each order of multipole- 
multipole interaction decays with differing power law behaviour. The agreement with the 
numerical Yukawa segment calculation is again good, except in the region around closest 
approach, where the platelets are co-axial, and the greatest deviation is observed, as seen 
in Fig.[|. This deviation decreases significantly the greater the vertical separation of the 
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platelets, and may be improved at close separation by the inclusion of higher order multipole 
moments. 

In Fig. |9], the angular dependence of the potential is examined by varying the angle 9 b 
at fixed centre-to-centre separation of r = 1.5a and 9a = 0. The charge-charge interaction 
is obviously independent of relative orientation, but significant variation in all higher order 
interactions is observed. Notably the quadrupole-quadrupole interaction is favoured when 
the discs are perpendicular to each other. The agreement with the computationally expensive 
numerical integration method is observed to be good, except when the discs are parallel, 
where the deviation is most noticeable, as observed and commented upon in Figs. 5 and 8. 

Finally in Fig. [10] the dependence of the interaction on the azimuthal angles is probed, 
for disc separations in the range 1 < r/a < 1.5 for discs fixed at 9a = 9b = 7r/4. As noted 
in Appendix B the interaction energy only depends on the difference A0 = (f>A — <Pb- As 
the separation of the discs increases the angular dependence of the pair potential is seen to 
diminish, and indeed at large separations disappears. 



VII. CONCLUSION 

The familiar multipole moment expansion of the electrostatic interaction between two 
extended charge distributions have been generalized to the case where the interaction is 
linearly screened by co- and counterions of an ionic solution in which the charged colloidal 
particles are dispersed. The case of uniformly charged, disc-like platelets, a model for the 
synthetic clay Laponite, was specifically considered in this paper, but the Yukawa-segment 
model, and the corresponding multipolar expansion, may be extended to particles of any 
shape. For example this procedure may readily be adapted to the case of uniformly charged 
rods, where symmetry again precludes multipole moments of odd order, by insertion of the 
relevant quadrupole and hexadecapole tensors for such a charge distribution. In the simplest 
case of spherical particles, the present treatment leads back to the familiar DLVO potential. 
For discs, the expansion, including up to quadrupole-quadrupole and charge-hexadecapole 
terms, yields interaction energies in good agreement with data for a discretized version of 
the Yukawa segment model, down to centre-to-centre distances of the order of the disc radius 
a, for all relative orientations of the two platelets which were investigated. As expected, the 
expansion breaks down at shorter distances, and yields rapidly divergent energies as r —>■ 0. 

The effective pair potential defined by eq.(pl|), and the explicit expressions in the Ap- 
pendices, should prove useful in Statistical Mechanics descriptions of semi-dilute clay disper- 
sions, and of their sol-gel transition, provided suitable short-range cut-offs are imposed. For 
computer simulations of more concentrated dispersions, an appropriate strategy would be 
to use a hybrid pair potential approach, interpolating between the multipolar expansion at 
large centre-to-centre distances, and a direct summation of the v 2 screened Coulomb site-site 
interactions in a discretized version of the Yukawa segment model, similar to that used in 



Ref. at short distances. 

In order to determine the phase behaviour of dispersions of charged platelets, from direct 
calculations of the free energy of systems of platelets interacting via the multipolar effective 
pair potential derived in this paper, it is important to include a structure-independent 
'volume' term in the free energy |fl8|. Such a volume term has been shown to play a crucial 



role in the determination of phase diagrams of suspensions of spherical charged colloidal 
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particles, in the regime of very low concentration of added electrolyte |||19|]. The exact form 
of the 'volume' term can be determined from a careful analysis of the density functional 
formulation of linearized PB theory [[§[15|. Such an analysis, which also provides a rigorous 
foundation of the Yukawa segment model fLl|, is under way. 
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APPENDIX A: CARTESIAN TENSORS FOR A SCREENED COULOMB 

INTERACTION 

The Cartesian tensors for a screened Coulomb interaction, are defined (analogously to 
the bare Coulomb case) as derivatives of the potential, ie: 

^... = V«V,V 7 ...(^^) (Al) 

Now, for any function f(r), (f(r) being e~ KT jr in the current work) the gradient V a /(r) = 
/ (j)V a r. For the bare Coulomb potential, f(r) = 1/r, and thus V a (l/r) = — (l/r 2 )V a r = 
T°. Combining these two simple results the gradient of a general function f(r) may be 
expressed in terms of the bare Coulomb tensor via V Q /(r) = — r 2 f {r)T®. Introducing the 
differential operator £>, defined by T>f(r) = —r 2 f (r), this may furthermore be written as 
V a /(r) = Vf{r)T®. Thus the set of successive interaction tensors for a screened Coulomb 
Potential T K may be written, using the product rule of differentiation as: 

K = V„/(r) = (p7)t q ° 

Tfr = V a V^(r) = (P 2 /)T Q %° + {V l f)T Q ap 

T K aPl = (V 3 f)T° Q T°T° + (V 2 f)[T° af3 T° + T a ° 7 T° + T^T° a ] + (V'fK^ 
T^s = (V 4 f)T° a T°T°T° 

+ (V 3 f)[T° p T°T° + Tppf + T° a5 T°T° + T» S T»T» + T° 7 T°T° + I* 2*Zj] 

_j_ (T)2 f\ rT~>0 rjiO I rpQ rpO i rpO rpO i rpQ rpO i rpQ rpQ . rpQ rpO i rpQ rpQ 1 

' \ u J )[ I aj3-y I 8 ' i a/35 i 7 + i a 7 <5 i /3 "+" 1 fl-yS 1 a + i a / 3 i 7 5 + I a-y I l3S "+" 1 a5 1 ^\ 

+ (V'fK^, (A2) 

where the coefficients {V n f} are functions solely of distance, given by 

V x j = (1 + Kr)e- Kr 

V 2 f = K 2 r z e~ Kr 

V 3 f = K 2 r\nr - 3)e- Kr 

V 4 f = K 2 r 5 (12 - 8/cr + K 2 r 2 )e~ Kr , (A3) 

and all factors of l/47reoe and more significantly all angular dependencies of the interaction 
are embodied in the bare Coulomb tensors T^a, expressions for which are easily calculated. 
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APPENDIX B: EVALUATION OF INTERACTION ENERGIES 



The interaction energy between two discs (A and B) , separated by a distance r AB and 
oriented at spherical polar angles (9 A ,(p A ) an d {9b,<Pb) may be written, eq. (|21"D , as a sum 
of contributions to the total from the interactions of each order of multipole via 



V AB (r A B, e A , 6 B , <Pa, <Pb) = Vfi + {Yab + O + Vg? + (Vjg + V%) + ■ ■ ■ , (Bl) 

where Z, Q and $ denote respectively the charge, quadrupole moment and hexadecapole 
moment on a disc, and V®b f° r instance denotes the contribution from the interaction of 
the charge on disc A with the quadrupole on disc B. 

The leading term in this expansion is simply the screened Coulomb interaction of the 
two charges, given by 



yff = Z A eT K Z B e 



Z 2 e 2 e~ Kr 
47ie e r 



(B2) 



where T K is the zeroth-order interaction tensor. The next term in eq. ( pi|) corresponds to 
the charge-quadrupole interaction, V^b , which is written in the following form 



v zq 

V AB 



-z a t:„q 



Ma/3 



(B3) 



2! 

Expressing the screened interaction tensor T£n in terms of the unscreened tensors {T®g...}, 
using eq.(|^2|), and recalling the definition of the quadrupole moment tensor, eq.(|13D this 
sum is calculated via: 



= ~Q [D 2 f[TX - T>f 2jnf] + V'f[T° a - nf 2>J ]} 



Q 



47reoe 
Za 2 1 



V 2 f 



1 — cos 2 9b 



+ V l f 



l-3cos 2 e B 



(B4) 



4 Ane^e r 3 



1 + nr + k r ) — cos 6 B {2> + 3nr + k r 



where {n B } denote the Cartesian components of the unit vector which defines the major axis 
of disc B, and the {£>"/} are defined by eq . (|A3|) . Along with this energy, the contribution 
due to the interaction between quadrupole on disc A and charge on disc B must be added. 
When all multiplicative factors have been included the final result reads 



V ZQ + QZ = yZQ+QZfo q Aj q b) 
Z 2 e 2 a 2 



47re e 
n 2 r 2 



x 



(B5) 



1 + /«r + — - (3 cos 2 6 A + 3 cos 2 B ) - 2 (l + Kr + kV) 



where it is observed that the charge-quadrupole interaction has no dependence on the az- 
imuthal angles 4> A and 4> B . 
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In the quadrupole-quadrupole interaction energy it is necessary to calculate the sum 
QapTap-ysQ^s- The screened tensor is expressed as a sum of terms involving the simpler 
unscreened tensors using eq. ( |A2[) , and simplified further by expressing all second and higher 
rank unscreened tensors in terms of the first order unscreened tensors, via: 



4ne eT^ = 15rT°T°T° - ^(T%, + Tp a , + T°5 Q/3 ). (B6) 

Following this procedure of expressing the elements of the n^-rank unscreened tensors in 
terms of those of the 1 st rank tensors, the full unscreened 4 th rank tensor, using ( |A2| ) and(|B6j), 
is given by 

T^s = T°T°T°T°[V 4 f + 18rV 3 f + 87r 2 V 2 f] 

+ [T°T°5 7 5 + T®T®5ps + T^Tgdp^ + T®T®5 a s + T^TgS ai + T®T$8 a p] 
[-(V 3 f/r 3 )-9(V 2 f/r 2 )] 

+ [<$ Q /3<$7<5 + SarySps + 5«$<W] [V 2 f/r 6 ] 

+ Vrfky, (B7) 



where the tensor T®p s appearing in the last line of eq. (|B7|) is the only term surviving if 
k = 0, corresponding to the purely Coulombic interaction. 

The screened Coulombic tensor T£p s , given by eq. (p7|) , is used to calculate both the 
contribution to the potential away from a single disc due to the hexadecapole moment 
^*(r) and also the quadrupole-quadrupole and charge- hexadecapole energies in the effective 
pair potential. For illustrative purposes the quadrupole-quadrupole interaction energy shall 
be pursued here, which involves the sum Qap^SpysQ^s- Using eq. (|E7|) it is evident that 
this sum will itself be the sum of contributions from terms involving Q A pT®T®T®T®Q^ 5 , 
QapT®Tp8 y sQjs etc. which must each be calculated separately. The first of these is calculated 
as 



Q*pT° Q T»pT»T»Q« s = Q 2 [8 a p - n A n A ][T»T»T»T»][8, s - n»nf] 

= Q 2 [T° a T° a - T»n A a T«n A ][T«T« - T°n»T°nf] 
Q 2 (l-cos 2 ^)(l-cos 2 ^) 
47re e r 8 

Proceeding along these lines for each of the terms appearing in QapT£p~sQ~s the interaction 
energy finally reads: 

" AB — 1^2^ap- L aP^5^i'yS 



Z 2 e 2 a A 

-(8 a p - n A n A )T^ s (5 lS - n*nf) 



64 

nsr (s^b) { [p4/ + 18rV3f + 87r2p2/] ^ (1 - cos2 6a){1 - cos2 6b) 



ii 



-[V 3 f + 9r£> 2 /] — [8 - 6 cos 2 9 A - 6 cos 2 9 B 

+4 cos 6U cos 9b (cos 9a cos #b + sin 9a sin #b cos(A0))] 
H — [6 + 2 (cos 9a cos #b + sin #a sin # b cos ( A0) ) 2 ] 

+ ^p[3 - 15 cos 2 9 A - 15 cos 2 9 B - 45 cos 2 9 A cos 2 # B 

+6(4 cos 9a cos #b — sin 6* a sin ^cos(A0)) 2 ]} (B9) 

where attention may be drawn to the fact that the interaction energy only involves the 
difference in the azimuthal angles, A0 = <p a — <t> b and not on their absolute values. 
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APPENDIX: TABLE CAPTIONS 

1. Coefficients appearing in the series expansion of the electrostatic potential along the 
z— axis, eq.(p|). C\ is the coefficient of the term of order x n in the binomial expansion 

2 

of (1 + x)* . 

APPENDIX: FIGURE CAPTIONS 

1. Integration over the surface of a disc 

2. Electrostatic potential along the z-axis. Solid lines denote the multipolar expansion, 
dashed lines the linearized PB potential, and triangles a numerical integration over 
a discretized charge distribution (Yukawa segment model). The upper set of curves 
are for na = 0.5, and the lower for Ka = 1.0. The divergence is highlighted in the 
logarithmic inset figure. 

3. Electrostatic potential in the xy-pleaie, symbols as in fig.(H) for na = 0.5,1.0. 

4. Geometry of a pair of platelets. The disc orientations are characterized by the spherical 
polar angles 9 and 0. The azimuthal angle 4>b has been omitted for clarity. 

5. Contributions to the Potential Energy as a function of separation for two parallel 
coaxial plates, (9a = 0,9b — 0, A0 = 0), for inverse screening length k = 0.5/a. The 
points plotted on the LPB curve (known in this geometry) correspond to a discretized 
Yukawa Segment model calculation. 

6. Contributions to the Potential Energy as a function of separation for two parallel 
coplanar plates, (9a = 7r/2,6>£ = 7r/2, Acf) = 0), again for na = 0.5. The discretized 
Yukawa segment calculation has been included for comparison. 

7. Contributions to the Potential Energy as a function of separation for two plates in 
T-shaped configuration, (9a = 0,9b = 7r/2,A0 = 0) with na = 0.5. The discretized 
Yukawa segment calculation has been included for comparison. 

8. Contributions to the Potential Energy as a function of the horizontal separation, r/a, 
for parallel plates with plate B fixed at a altitude of 1.5a above disc A and na = 0.5. 
The discretized Yukawa segment calculation has been included for comparison. 

9. Contributions to the Potential Energy as a function of the angle 9b with 9a = 
for na = 0.5, and fixed centre-to-centre distance r = 1.5a. The discretized Yukawa 
segment calculation has been included for comparison. 

10. Dependence of the energy on the difference in azimuthal angle A</>, plotted for fixed 
9a = 9b = vr/4 for disc separations r/a = 1.0,1.1,1.2,1.3,1.4 and 1.5 from top to 
bottom, for na = 0.5. 
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n 









+1 





1 
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2 






4 


' 24 


3 





1 K« 4 K 3 a 6 
8 192 


A 


+ 8 


K 2 a 6 , K 4 a 8 
32 1 1920 


5 







64 1 192 23040 


n even 


2a n C (n/2) + l 

2 




n odd 
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